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AROUND A BLACK HOLE 
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Abstract. We solve a class of boundary value problems for the stationary ax- 
isymmetric Einstein equations corresponding to a disk of dust rotating uniformly 
around a central black hole. The solutions are given explicitly in terms of theta 
functions on a family of hyperelliptic Riemann surfaces of genus 4. In the absence 
of a disk, they reduce to the Kerr black hole. In the absence of a black hole, they 
reduce to the Neugebauer-Meinel disk. 
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1. Introduction 



Two of the most famous solutions of the stationary axisymmetric Einstein equa- 
tions are the Kerr black hole and the Neugebauer-Meinel disk. The former was 
discovered by Kerr in 1963 |Kej and the latter by Neugebauer and Meinel in the 
1990s [NMH INM21 INM3j . In this paper, we construct analytic solutions of a class 
of BVPs for the stationary axisymmetric Einstein equations which combine the Kerr 
and Neugebauer-Meinel spacetimes. Thus, the BVPs considered involve a finite disk 
of dust rotating uniformly around a central black hole. In the limit of a vanishing 
disk, the solutions tend to the Kerr black hole. In the absence of a black hole, they 
reduce to the Neugebauer-Meinel disk. The constructed disk/black- hole systems are 
given explicitly in terms of theta functions on a family of hyperelliptic Riemann sur- 
faces of genus 4. Given the importance of the Kerr and Neugebauer-Meinel solutions, 
we believe that the class of solutions presented here could also be of interest. 

The general analysis of rotating relativistic bodies is exceedingly complicated be- 
cause it involves the study of free boundary value problems (BVPs) for the Einstein 
equations, which are nonlinear partial differential equations in four dimensions. How- 
ever, in cases where the surface of the body is known and the motion is stationary and 
axisymmetric (a reasonable assumption in many astrophysical situations), the physi- 
cal problem gives rise to a BVP for a single integrable equation in two dimensions — the 
celebrated Ernst equation. The integrability of the Ernst equation implies that pow- 
erful solution-generating techniques are at hand. Thus, through the application of 
suitable nonlinear transformations, new stationary axisymmetric spacetimes can be 
generated from already known ones. Furthermore, a large class of solutions of the 
Ernst equation can be given explicitly in terms of theta functions on Riemann sur- 
faces 5K|. In this way, it is possible to write down a large number of exact analytic 
solutions to the stationary axisymmetric Einstein equations and to study them using 
the methods of algebraic geometry cf. |KR2j . 

Nevertheless, for the solution of a concrete BVP, the power of this approach is often 
limited. Indeed, although a large class of exact solutions can be produced, the problem 
of determining which particular solution within this class that satisfies the given BVP 
is in general a highly nonlinear problem. It is therefore remarkable that Neugebauer 
and Meinel in the 1990s were able to solve explicitly, using constructive methods, 
the BVP corresponding to the physically relevant situation of a rotating dust disk. 
The structure of an infinitesimally thin, relativistic disk of dust particles which rotate 
uniformly around a common center was first explored numerically by Bardeen and 
Wagoner 40 years ago [ BW2j . who also pointed out that "there may be some hope of 
finding an analytic solution" [BW| . Bicak notes in the comprehensive review pBj that 
the subsequent construction of such an analytic solution by Neugebauer and Meinel 
represents "the first example of solving the BVP for a rotating object in Einstein's 
theory by analytic methods." Let us point out that since the Neugebauer-Meinel 
dust disk can be written in terms of theta functions on Riemann surfaces of genus 
2, it belongs to the general class of solutions introduced in [K] and can therefore be 
analyzed by means of algebro-geometric methods. On the other hand, the Kerr black 
hole is the most famous example of a stationary axisymmetric spacetime and has had 
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an immense impact on the development of general relativity and astrophysics (see 
e.g. P). 

The approach in this paper is primarily inspired by the work of Neugebauer, Meinel, 
and collaborators |MAKj . but also by a novel method for the analysis of BVPs for 
integrable PDEs which has been developed by Fokas and his collaborators within the 
past 15 years. A central development in the theory of nonlinear PDEs in the second 
half of the 20th century, and continuing to the present, has been the introduction 
of the Inverse Scattering Transform (1ST). This technique was put forward in the 
famous 1967 paper |GGKj in connection with the Korteweg-de Vries (KdV) equation 
and the range of its applicability began to unfold with the investigation of the nonlin- 
ear Schrodinger (NLS) equation [ZS]. One of the most important later developments 
in this area has been the generalization of the 1ST formalism from initial-value to 
initial-boundary value problems introduced by Fokas \FT\ IF2j and subsequently de- 
veloped further by several authors cf. |F3j . The Fokas method consists of two steps: 

(a) Construct an integral representation of the solution characterized via a matrix 
Riemann-Hilbert (RH) problem formulated in the complex /c-plane, where k denotes 
the spectral parameter of the associated Lax pair. Since this representation involves, 
in general, some unknown boundary values, the solution formula is not yet effective. 

(b) Characterize the unknown boundary values by analyzing a certain equation called 
the global relation. In general, this characterization involves the solution of a nonlin- 
ear problem; however, for certain so-called linearizable boundary conditions, step (b) 
can be solved in closed form. 

In a recent work [LFj steps (a) and (b) were implemented for the class of BVPs of 
the Ernst equation corresponding to a thin rotating disk of finite radius. In particular, 
it was found that the boundary conditions of the particular EVP corresponding to 
the uniformly rotating Neugebauer-Meinel disk are linearizable. The present paper is 
to some extent a continuation of [LFj : The main novel observation is that the EVP 
which combines the Kerr black hole boundary conditions with those of a uniformly 
rotating disk is also linearizable. 

Physically, disk/black-hole systems are important as models for black hole accretion 
disks in the context of active galactic nuclei and X-ray binaries [A]. Accretion disks 
are flattened astronomical objects made of rapidly rotating gas which slowly spirals 
onto a central gravitating body. Accretion onto a black hole is generally assumed to 
be thin and axisymmetric (see e.g. [Pj) and many of the most energetic phenomena 
in the universe have been attributed to the accretion of matter onto black holes. 
In particular, active galactic nuclei and quasars are believed to be the accretion 
disks of supermassive black holes. We refer to and references therein for more 
information on the physical aspects of black hole accretion disks. Our approach here 
is mathematical and we do not investigate any possible physical relevance of the 
presented solutions. 

The manuscript is organized as follows. In section [2] we formulate our disk/black- 
hole EVP and write down its full solution in terms of theta functions. In section [3j 
we consider a particular example. In sections 4][7, the derivation of the solution is 
presented in several steps. In section [8j we consider the singularity structure of the 
solution and its dependence on various parameters. In the appendix, we consider the 
relationship between the solution derived here and the general class of solutions of 
the Ernst equation studied in [KMl IKKSj . 
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2. Disk/black-hole systems 

In this section we introduce the Ernst equation, formulate the BVP corresponding 
to a dust disk rotating uniformly around a central black hole, and present its solution 
in terms of theta functions. 

2.1. The Ernst equation. The metric of a stationary axisymmetric vacuum space- 
time can be written in the Weyl-Lewis-Papapetrou form 

ds^ = [e^'^idp^ + dC^) + p^dif^] - e^^'idt + ad^pf, (2.1) 

where p >0 and C G IK are Weyl's canonical coordinates and t, ip are chosen so that dt 
and dip are the two commuting asymptotically timelike and spacelike Killing vectors, 
respectively, cf. |SKMj . The metric functions e^^ , a, and e^'^ are functions of p and 
C alone. Introducing a real-valued potential b by 

ap = pe-^^h^, = -pe-^^bp, (2.2) 



it can be shown that the Einstein field equations for the metric (2.1) reduce to the 
following single nonlinear PDE in two dimensions for the complex-valued Ernst po- 
tential f = e^^ + ib: 



f + f 
2 



/pp + /« + ^/p) = /p + /O P>0, CGK. (2.3) 



Following standard practice, the real part of / will be denoted by e'^^ despite the fact 
that it may take on negative values. 

We also need the concept of a corotating frame. Given G M, we define the 
coordinates {p' ,C,' ,(p\t') corotating with the angular velocity Q by 

p = p^ C' = C) = V ~ ^t, t' = t. 



In these new coordinates, the metric (2.1) retains its form and the corotating metric 



functions UQ,aQ, kq are related to U, a, k via 

^2Un ^ ^2U ^ f^^)2 _ ^2^2^-4U^ ^ (2.4a) 

(1 - nan)e^^'' = {1 + na)e^^ , nn-Un = K-U. (2.4b) 

The Ernst equation retains its form in the corotating system and we denote the 
corotating Ernst potential hy fn := e"^^^ + ibn. 

2.2. The boundary value problem. We now formulate the BVP for the Ernst 



equation (2.3) which corresponds to a finite dust disk rotating uniformly around a 
central black hole. The formulation involves five parameters: the radius po > and 
angular velocity f2 G M of the disk; the 'radius' ri > and angular velocity O/j G M of 
the black hole horizon; and the (necessarily constant) value of the corotating metric 
function e^^" on the disk. We will henceforth work with a complex variable z = p+iC, 
and write f{z) for /(/O, C)Q 

Let P denote the exterior of a finite disk of radius po > 0, i.e. T) consists of all 
z G C with strictly positive real part which do not belong to the interval [0,/0o], see 

^In general, given a function h{p,(^), we will suppress the dependence on z and write h{z) for 
h{p, (^) even when h is not analytic. 
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Figure 2.1. The exterior domain D of a disk of radius po with the black hole horizon 
stretching along the imaginary axis from —iri to iri . 



Figure 2.1 The horizon of the black hole stretches in the z-plane along the imaginary 
axis from —iri to iri. We consider the problem of finding a function / such that: 



/ satisfies (2.3) in T? 



f{z) = f{z) (equatorial symmetry). 

f{z) — )■ 1 as — )■ oo (asymptotic flatness). 



(2.6a) 
(2.6b) 
(2.6c) 



dp 



{iC) = for all \C, \ > ri (regularity on the rotation axis) 



f'Q,{p ± iO) = e^^"(+*°) for < p < po (boundary condition on the disk). (2.6d) 
g2t/f2^(«C) _ Q £qj. < |(^| < ri (boundary condition on the horizon). (2.6e) 



The boundary conditions (2.6d) and ( 2.6e| ) are the boundary conditions corresponding 
physically to a uniformly rotating dust disk and a rotating black hole, respectively, 
cf. |MAK| . If one sets ri = in (2.6) (i.e. one removes the black hole), then the 
solution of the obtained BVP is the Neugebauer-Meinel disk rotating with angular 
velocity fi. If one sets po = in (2.6) (i.e. one removes the disk), then the solution 
of the obtained BVP is the Kerr black hole rotating with angular velocity Qh- 



2.3. The solution. The formulation of the BVP (2.6) involves the five independent 
parameters po, f^, ri, VLh-, and the constant value of e^^" on the disk. However, it turns 
out that the condition that the solution be nonsingular at the rim of the disk imposes 
one relation among these parameters, so that the class of solutions is parametrized 
by only four parameters. It is convenient to adopt a parametrization in terms of 
the four parameters po, ri, W2, and w^, where W2 and w/i are two real quantities 



related to the other parameters via equations (5.9) below. For a given choice of these 



parameters, the corresponding solution / of the BVP (2.6) can be written in terms 



of theta functions on the family of Riemann surfaces {S^j^gx) defined as follows: Let 
fci, ^1, . . . , ^4, ^4 G C denote the eight zeros of w'^{k) + 1, where 



w{k) 



w^k^ + W2k'^ + Poiw2 - W4pI) 
{k^ - rf) ' 



(2.7) 



ordered so that /cj, j = 1, . . . , 4, have negative imaginary parts and so that 



Re ki < Re ^2 < Re k^ < Re k^. 
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-IPo 



IZ 



-iz 



Figure 2.2. T/ie Riemann surface T,z presented as a two-sheeted cover of the complex 
k-plane with five branch cuts when Rek^ < C- The contours 7 and are also shown. 



Since these eight zeros are symmetrically distributed with respect to the origin, we 
have —ki = k^ and —k2 = ^3. For each z = p + iC,, T,z is defined as the hyperelliptic 
Riemann surface of genus 4 consisting of all points (/;;, y) G such that 

4 



= {k + iz){k - iz) W{k- kj){k 



(2. 



together with two points at infinity required to make the surface compact. We intro- 
duce branch cuts in the complex /c-plane from kj to kj, j = 1, ... ,4, and from —iz 
to iz, see Figure 2.2 For fc G C = C U {00}, we let k^ and k~ denote the points 
which project onto k and which lie in the upper and lower sheet of E^, respectively. 
By definition, the upper (lower) sheet is characterized by y/k^ — )• 1 {y/k^ — )• —1) as 
/c — )• 00. 



In view of the assumption of equatorial symmetry (2.6b), we will in the sequel 



always assume that C ^ 0- We will also, for the sake of definiteness, assume that 

< Re/c3 < n < ReA;4. (2.9) 

Moreover, we assume that ^ 7^ Re/ca and C ^^k^, so that no branch cuts overlap 
(the solution for ( = Kekj, j = 3,4, can be obtained by continuity). 



For n complex numbers {aj}", we let [ai, 



denote the directed contour 



1 n 1 



Oj+i]. We let 7 denote the contour on which projects to the contour 



[n , Re ks + e, fcg + e] U [ks - e. Re ks - e, Re /c2 + e,k2 + e]U [k2 - e, Re k2 - e, -n] (2.10) 

in the complex fc-plane, where e > is an infinitesimally small positive number, and 
which lies in the upper sheet for Re A; < ^ ^md in the lower sheet for Re A; > C. 
We define F^ as the contour in the upper sheet of which lies above the segment 
F = [—ipQjipo]. The contours F+ and 7 are shown in Figure 2.2 in the case when 
Refc4 < C- 

In order to define theta functions associated with S^, we need to introduce a basis 
of the first homology group i7i(S2,Z) of S^. Since the Riemann surface depends 
on z, there are three qualitatively different cases to consider: 

(1) The cut [—iz,iz] lies to the right of [k^, k/i] (i.e. Rek^ < (). 

(2) The cut [—iz,iz] lies between [^3,^:3] and [A;4,fc4] (i.e. Re/c3 < C < ReA;4). 

(3) The cut [—iz,iz] lies to the left of [/i:3,A;3] (i.e. < C < ReA;3). 
For these three cases, we let {aj, ^j}'^=i be the canonical basis of Hi(T,z,1j) shown in 

1, 



1), (2), and (3) of Figure 2.3, respectively. Thus, for j 



,4, Oj surrounds the 



BOUNDARY VALUE PROBLEMS FOR THE EINSTEIN EQUATIONS 7 




of genus g = A in the case of (1) Rek4 < (2) Reks < ( < ReA;4, and (3) 
< C < ReA;3. 




Figure 2.4. Three-dimensional picture of the hyperelliptic Riemann surface T^z in 
the case when Rek^ < i^. 
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cut [kj,kj], whereas bj enters the upper sheet on the right side of [—iz,iz] and exits 
again on the right side of [kj, kj]. 

We define {wj}^ as the canonical basis of the space of holomorphic one- forms on 
dual to {aj,bj}. Then 

u}i = 6ij, / uji = Bij, i,j = l,...,4, 

■J ''^3 

where B is the period matrix associated with the cut system {aj,6j}. We let uj = 
(wi, 072, W3, a;4)'^. The 4x4 matrix B is symmetric and has a positively definite 
imaginary part. The associated theta function Q{v\B) is defined by 

Q{v\B) = e2-(l^^^^+^^''), V G (2.11) 

We let ujpQ denote the Abelian differential of the third kind on S^, which has two 
simple poles at the points P and Q with residues +1 and —1, respectively, and whose 
a-periods vanish, i.e. ^pq = for j = 1, . . . , 4. 
We can now state our main result. 

Theorem 2.5 (Solution of the disk/black-hole BVP). Let po,'i^i,W2,W4 be strictly 
positive numbers such that (S.Sl^ holds. [The requirement that the solution be singularity- 
free imposes further restrictions on these parameters, see section^] Let the function 
h{k) be defined by 

h(k) = ^ In [^Jw{kf + 1 - w{k)) , keV =[-ipQ,ipQ]. (2.12) 
Define the z-dependent quantities n G and I £M. by 



hoj+ I UJ, I = J^^hu}^+^- + J t^oo+oo-- (2.13) 



Then the function 



/(^) = -7 (2.14) 



satisfies the BVP (2.6) with the prescribed values of pQ and ri, and with the values 
ofQh, ^, and 6^*^"^+*^ given by 



(2.15) 

where Uhor £ 1^ denotes the (necessarily constant) value of the metric function a on 
the horizon and e^^° G M denotes the real part of f{-\-iO). Explicit expressions for the 
constants Uhor and f{+iO) are presented in propositions 2.9 and 2.10 below. 
Moreover, define the z-dependent quantity L by 



hoo + - 
r+ "i"! 

(2.16) 



^ lim ( / UJ + - - sgn(C - n) [ uj + - 
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where 71(e) denotes the contour 7 with the segment covering [—ri, —ri + e] removed, 
72(e) denotes the contour 7 with the segment covering [ri — e, ri] removed, and the 
prime on the integral along T indicates that the integration contour should be deformed 
slightly before evaluation so that the pole at K 2 = is avoided^ Then the metric 



functions e , a, e of the line element (2.1) corresponding to the Ernst potential 



(2.14) are given for z £1) by 



0(0) I 

Q{u) ' 



a[z] 



ao 



Q(0) 



g(o)e(u + /ia;|i?) 



Q{u) 



eiu\B)Qiu + rJ^^uj\B) 



e{0\B)Q{Zlu;\B) 
where Q{v) is defined by 



Q{v) 



uj\B) 



e[v\B)Q{v + PX^\B) 



V e 



(2.17a) 
(2.17b) 

(2.18) 



and the two constants ao,Ko G 
below. 



are given explicitly by equations (2.26) and (2.21) 



Remark 2.6. 



1. The function h{k), A; G F, is well-defined by the right-hand side of 

w > Q for A; G r. 



(2.12), because + 1>Q and ^/w 



+ 1 



2. Unless stated otherwise, all integrals in this paper along paths on Riemann 
surfaces for which only the endpoints are specified are assumed to lie within the 
fundamental polygon obtained by cutting the Riemann surface along the given cut 
basis. This implies that some integrals will depend on the particular choice of the 
a,'s and 6j's within their respective homology classes. It is convenient to fix the Oj's 



and 6j's so that they are invariant under the involution 
is accomplished by fixing^' 



k^. For C > Re A;4, this 



64 = [iz, k/i\^ U [/c4, iz] 
and by letting Cj, j 



•J - bj+i U [kj+i,kj] ■ u [Kj,Kj+n 

4, be the path in the homology class specified by Figure 

which as a point set consists of the points of lying directly above [kj, kj]. This 

implies the important symmetry u}{k~^) = —u}{k~). For with < (" < Re/c4 and 

for other Riemann surfaces below, we will assume that an analogous fixing of the cut 

basis which assures u{k~^) = —uj{k~) has been made. 

3. The limit as e — )• of the expression within brackets on the right-hand side of 

and uj,.+„- . 
1 '1 



U[kj,kj+i] , j = 1,2,3, 



1, 



2.3 



(2.16) always exists and is finite because of the pole structure of a;_^+ 

1 



4. The assumption that ri satisfies (2.9) is made for simplicity and is not essential. 



The relevant formulas in the case when (2.9) does not hold can be obtained from 



those presented here by analytic continuation. 

2.4. Axis and horizon values. Of particular interest are the values of the Ernst 
potential and of the metric functions on the regular axis {i^ \C > '^i} ^ind on the 
black hole horizon {^C I < C < ^i}- Iii the limit /) — )• 0, the Riemann surface 
degenerates since the branch cut [—iz, iz] shrinks to a point. Thus, the values of / on 



The result is indepedent of whether the contour is deformed to the right or to the left of the pole. 



'^For two complex numbers zi and Z2, [zi,Z2j 
upper and lower sheets of , respectively. 



and 



21,^2 



denote the covers of [zi , 22] in the 
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the C-axis are given in terms of quantities defined on the Riemann surface S' defined 
by the equation 

4 

y'^ = Y[(k-k,){k-k,), (2.19) 

i=i 

i.e., S' is the Riemann surface with the cut \—iz, iz] r emoved. We introduce a 



canonical cut basis {a^-,6^}f on S' according to Figure 2.7 and let co' = {ljj[,ljJ2,uJ' 



denote the dual basis of holomorphic one- forms. We let B' denote the associated 
period matrix and introduce the short-hand notation &'{v) := @{v\B'), v G C^, for 
the associated theta function. 



Let 7' denote the C-dependent contour on S' which projects to the contour (2.10) 
in the complex A:-plane and which lies in the upper sheet for Re k < ( and in the 
lower sheet for Refc > C- Define the ("-dependent quantities u' G C^, /' G M, and 

K' eChy 

u = [ huj' + [ J, I' = [ hJ^+^.+ j = [ ^c+C-- 

(2.20) 

Let 7"'' denote the contour on S' with the same projection onto the complex /c-plane 
as 7', but which lies entirely in the upper sheet, and define J' G M by 

The prime on the integral along 7"^ indicates that the path should be deformed 
slightly before evaluation, so that it avoids the pole of the integrand at A; = C^- It 
is irrelevant for the formulas below if this deformation is performed so that the pole 
lies to the right or to the left of 7"^, since these two choices yield values of J' which 
differ by a multiple of 27rz and J' only appears exponentiated. In the same way, the 
value of J' changes by irrelevant multiples of liri if loops surrounding C,^ are added 
to the contour from rj~ to . 



Proposition 2.8 (Solution on the regular axis). The behavior of the solution (2.14) 
near the regular axis {iC | C > ''"'1} given by 

f{p + iC) = f{iC) + 0{p^), p^O, C>ri, (2.22) 
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where 



m) = . 1- . . 1- C>n. (2.23) 



The behavior of the metric functions e^^ , a, and e'^'^ in (2.11) near the regular axis 
is given by 

g2l/(p+^C) = ^2U{^0 + o{p^), a{p + iC) = 0{p^), 



p^o, C>n, 



where 







J'-J 


e'(o)2 







C > n- 

(2.24) 



Define the constant Lq by 

Lo = -\j^dn,'^^{K,) j^h{n,)u:'^,^.{^^^ (2.25) 



H — lim 
2 e^O 



OJ 



2 In e , 



72^ {^) 



where denotes the contour 7"*" with the segment covering [— ri, — ri+e] removed, 

and 72" (e) denotes the contour 7''' with the segment covering [ri — e,ri] removed. 
Define the constants uq^Kq G M by 



qq = — 2i 



Q'iu' 



e ^ lim ( i?e ^4 00+00- ) 



(2.26) 
(2.27) 



where the right-hand sides are understood to be evaluated at some C > ''^10 Define 
the ("-dependent quantities L' and M' by 

L' = -- [ dKi^(Ki) [ h(K2)uj'+ -(k+)+ / huj' + _ + / W+ _ (2.28) 

< C < r-i, 



H — lim 

2 e^O 



and 



M' = - lim 

2 e^O 



7((e) '1 



21ne 



r+ 



(C-^)- 



c+c- 



21ne-ln4-7ri , < ( < n, (2.29) 



where 7[(e) denotes the contour 7' with the segment covering [— ri, — ri + e] removed, 
and 72(e) denotes the contour 7' with the segment covering [ri — e,ri] removed. 



Proposition 2.9 (Solution on the horizon). The behavior of the solution {2.14) near 
the black hole horizon | < ( < ri} is given by 

f{p + K) = f{K) + 0{p^), p^O, 0<C<ri, (2.30) 



^Their values are independent of the choice of ^ > 
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where 



m) 



< c < ri. 

(2.31) 

The behavior of the metric functions e^^ , a, and e^'^ in (2.11) near the black hole 
horizon is given by 



2K(p+iC) _ 



+ 0(p2), a{p + i0 = ahor + O{p^), 
+ 0(p2), p^o, 0<C<ri, 



where 



G'(0)2 



(2.32) 



(2.33) 



J'+J' 



Oo + 



e'(u' + 2/°^"L^')0'(o)^ 



-I'-M' 



Q'{u' + f^^ 00') (O'if^' oj')^ - Q'if^' a;')2e-2^')^ 
' G'(0)2 

Moreover, Ohor O'f^d e^'*''"'' are constants independent of < ( < ri. 



, (2.34) 
(2.35) 



By taking limits in the above formulas we can find the values of / at the origin 
z = +iO and at the point z = iri where the regular axis meets the horizon. 

Proposition 2.10 (Solution at z = +iO and z = iri). The value of the solution 
(2. 14-) at z = iri is given by 



fiin) 



GV + J.-'c.')' 



(2.36) 



C=ri 



Define the value of J' at = by 



J'k=o 



+ 



0+0- ' 



where the primes on the integrals indicate that the contours should be deformed slightly 
before evaluation so that they pass to the left of the pole at k = 0'^ . Then the value 
f{+iO) of the solution (2.14) at z = +iO is given by the right-hand side of (2.31) 
evaluated at ( = 0. 



3. Example 

Before presenting the derivation of the solution presented in the previous section, 
we wish to consider a concrete example. In this regard we note that all the quantities 
appearing in section [2] can easily be computed explicitly using standard results for 
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Riemann surfaces. In order to find the canonical basis ^'^d the period matrix 

B, we start with the holomorphic one-forms {Oli defined by 

Q = -^. i = l,...,4. (3.1) 

y 

The Ci's form a noncanonical basis of holomorphic differentials (see e.g. jFKj ) . Defin- 
ing the two 4x4 matrices A and Z by 



= [ 







^,3 = 1, 



we find 

The one-form oo, 



u = AC, B = AZ. 
oo^oo- on is given exphcitly by 



whereas the one- form ^^+^- on S' is given by 



k^dk 



(3.2) 
(3.3) 

(3.4) 



y'{r)dk 

(k-Ov'ik) 



.[k-Q)y'{k) ^ 



(3.5) 



For /c's which lie some distance away from the branch cuts, the value of y in (2.8) can 
be evaluated according to 

4 



y{z, k+) = {k + tz)\j \[{k - k, 



Ik-k 



'3 



k-k 



k G 



'J 



where the branches with strictly positive real part are chosen for the square roots. 
For k £ [—iz, iz] and e > infinitesimally small, we have 



yiz,ik + e)+) = (k + iz) -i. 



k — iz 



k + iz 



Hik-k,) 



k kj 



k E 



similar expressions are valid when k S [kj,kj]. Using formulas of this type, it is 
straightforward to numerically evaluate all the expressions presented in section [2j In 
fact, the theta functions are particularly suitable for numerical evaluation, because 
the strictly positive imaginary part of B implies that only a small number of terms in 
the sum (2.11 ) have to be included. In this way, we have verified for several examples 
all formulas of section [2] to high precision. 

3.1. Numerical data. Here we consider the particular example for which 
Numerically, we find 

A;i w -0.95 - 5.48i, A;2 -0.21 - 0.95i. 



Po 



ri 



(3.6) 



Thus, (2.9) is satisfied and we compute 

n 0.055, ~ 0.14, f{+iO) 

and 

ao ~ -18.17, Ko ^ 9.43, 



-0.17, 



f{iri) s 
-93.46. 



-OMi, 
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The real and imaginary parts of / are shown in Figure 3.1 Note that the real part 



of / is continuous but not smooth at the endpoints ibiri of the horizon. Moreover, 
as expected, the imaginary part of / has a jump across the disk. The values of / in 

the equatorial plane = 0"*" are shown in Figure 3.2 

The metric functions ae^^ and e'^'^ are shown in Figure 3.4 Note that ae^^ 



on 



the regular part of the axis. Moreover, e = 1 on the regular axis and e 
on the horizon. Figure 3.5 shows the real parts of the corotating potentials fn and 
/q^ in the equatorial plane and along the C-axis, respectively. In accordance with 
the boundary conditions (2.6d) and (2.6e), e^^" is constant along the disk and e^^^h 
vanishes along the horizon. It can also be verified to high accuracy that the metric 
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Figure 3.4. The metric functions ae^^ and e^'^ for the choice of parameters specified 



in (3.6). 




0.05 r 



Figure 3.5. The real parts of the corotating potentials fn and f^^^ in the equatorial 
plane and along the (-axis, respectively. 



functions a and k defined by (2.17) satisfy tlie appropriate equations, i.e. (cf. |KR2| ) 



1^1 



and Kz 



2e^u 



fzfzj 



and that 



doQ 



dC 



< p < po- 



(3.7) 



(3.8) 



z=p+iO 



Equation (3.8) implies that the imaginary part of is constant along the disk in 
accordance with (2.6d). 



4. Spectral theory 

The now turn to the proof of the results of section [2j The proof will proceed 
through four main steps, presented in sections |4j [5j JHland [7| respectively. The first 
step consists of analyzing the Lax pair for equation ( |2.3[ ) and formulating two matrix 
Riemann-Hilbert (RH) problems: one main RH problem (which can be formulated 
for any choice of boundary conditions) and one auxiliary RH problem (which can be 



formulated because the boundary conditions of the BVP (2.6) are linearizable). In 
the second step, we show that these two matrix RH problems can be combined into a 
single scalar RH problem. The third step consists of solving this scalar RH problem 
explicitly in terms of theta functions. In the final fourth step, we prove propositions 
2.8||2.10 concerning the behavior of the solution near the (^-axis; this step will follow 



from a study of the theta function formulas of theorem 2.5 in the limit p — )• 0. 
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4.1. A bounded and analytic eigenfunction. The elliptic Ernst equation (2.3) 
admits the Lax pair 

U,{z,k) = U{z,kMz,k), 
\Mz,k) = V{z,k)^z,k), 

where the 2 x 2-matrix valued function ^{z, k) is an eigenfunction, is a spectral 
parameter, and the 2 x 2-matrix valued functions U and V are defined by 



We write the Lax pair (4.1) in the form 

= where W := Udz + Vdz. (4.2) 

Let 

l\ fO -i\ fl 



Suppose / is a solution of the BVP (2.6). Following the same procedure as in |LF| . 
we define a solution ^{z, k) of (4.2) with the following properties: 

• For each z, ^{z, •) is a map from the Riemann surface Sz to the space of 2 x 2 
matrices, where Sz is defined by the equation 

^ = 1 ^• 

k + iz 

We view 5^ as a two-sheeted covering of the Riemann A;-sphere endowed with 
a branch cut from —iz to iz; the upper (lower) sheet is characterized by A — )• 1 
(A — )> —1) as — )• oo. 

• $ satisfies the initial conditions 

hm mz,k')\i = fj") , lim [^{z,k+)]2 = f , A; G C, (4.3) 

where, for a 2 x 2-matrix A, [A\i and [A\2 denote the first and second columns 
of ^4, respectively. 

• <I> obeys the symmetries 



^>(z, A:+) = o-3$(z,A;")cri, ^{z,k+) = (Ti^{z,k+)a^, keC. (4.4) 

• <I> is analytic for k £ Sz away from the set U F^ U {—r^Tr^}, where F"^ 
and F~ denote the coverings of F = [—ipo,ipo] in the upper and lower sheets 
of respectively. 

We emphasize that in general, has singularities (simple poles) at the points k^ for 
k = — ri and k = ri. These poles arise since the Lax pair (4.1) is singular at points 
where e^^ = 0. Physically, the points at which e^^ vanishes make up the ergospheres 
of the spacetime (within these surfaces there can be no static observer with respect 
to infinity). To see that e^^ = at z = iri and z = —iri, we note that the metric 
function a vanishes identically on the regular axis (cf. |SKMj ). Thus, evaluating 
(2.4a) at p = 0, we find Ufi{iC) = U{iC,) for C, > n. The boundary condition (2.6e) 
together with the continuity of imply that e^^^**"^ — )■ as C i ^'i • 

In addition to the eigenfunction ^[z^k), we will also need its corotating analog 
$0(2:, /c) . This eigenfunction satisfies the Lax pair equations (4.1) with / replaced 
by fn and the initial conditions (4.3) with $ replaced by ^q,. The eigenfunctions ^ 
and are related by jMAK j 

^n{z,k) = Kn{z,k)^{z,k), k G Sz, (4.5) 
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where 

A^{z, k) = {l + na)I - npe-'^'^as + i{k + iz)ne-'^'^ {X{z, k)ai - 1)0-3 (4.6) 

and I denotes the 2x2 identity matrix. The corotating eigenfunction ^Q^^{z,k) is 
defined analogously. 

4.2. The main Riemann-Hilbert problem. Evaluation at /> = of equation 
(2.4a) with replaced by il/i yields 

e2^%«) = e2^(^^)(l + ^^/,a(^C))^ < C < n. (4.7) 

Note that e^^ < along the horizon, which lies inside the ergosphere. Thus, in view 
of the boundary condition (2.6e), we find that a = ahor on the black hole horizon, 
where ahor is a constant given by 

ahor = -l/^h- (4.8) 

The next proposition expresses the values of <^ on the ^-axis in terms of two spectral 
functions F{k) and G{k). We let /i denote the value of / at 2; = iri. 

Proposition 4.1. The values of <I> on the Q-axis can he expressed in terms of two 
spectral functions F{k) and G{k) as 



^«,^+)= ^^•^"^^ 

^«,^+)=(^[|j }^Ti{k)Aik), 0<C<ri, keC, (4.9b) 

HiC, k+) = (^j^^l _\) T2{k)aiAik)au -n < C < 0, keC, (4.9c) 

HiC,k+) = (^^l^^j }^ cJiA{k)ai, C < -n, keC, (4.9d) 
where the 2 x 2-matrix valued functions A{k), Ti{k), and T2{k) are defined by 

Ti{k) = ^r^^('^^''-'''^.%'-'^' Ao^-.V k€C, (4.11) 

^ 2{k-ri)nh\ 2{k-ri)nh + ifiJ ^ ' 

T2{k) =Tii-k), keC. 

The functions F{k) and G{k) have the following properties: 

• F and G are unique functions of k £ C, i.e. viewed as functions on Sz they 
satisfy 

F{k+) = F{k-), G{k+) = G{k-), keC. (4.12) 

• F{k) and G{k) are analytic for A; G C \ (F U {ri, — ri}). 

• Under the conjugation k ^k, F and G obey the symmetries 



F{k) = F{k), G{k) = -G{k), keC. (4.13) 
In the limit A; — )• oo, 

F{k) = l + Oil/k), G{k) = 0{l/k), k^oo. (4.14) 
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Proof. For z = X = 1 for all k on the upper sheet of Sz- The axis values of ^ 
are thus determined by integration of the equation 



d<^ = W{iC,k+)^ 



where 



W{iC,k+) 



df df 



f + f\df df)- 



(4.15) 



(4.16) 



Since the real part of / vanishes at 2; = itzri, equation (4.15 ) breaks down at (" = i^'i- 
Integration of (4.15) for Q in each of the four intervals (— 00, — ri), (— ri,0), (0,ri), 
and (ri, 00) yields 



^«,A:+) 
^iCk+) 
^«,A:+) 
$(iC,A:+) 





1 


/«) 




/«) 








fiK) 




fiiC) 













C>ri, /ce 

< C < n, G 

-ri < C < 0, A; e 

C < -ri, A; G 



where the matrices Uj{k), j = 1, ... ,4, are independent of C- The initial conditions 
(4.3) imply that Ui = A for some functions F{k) and G{k). This establishes (4.9a). 
The value of ^> at 2; = —zoo is obtained from the value at z = ioo by integrating 
along a large semicircle at infinity. During this integration k changes sheets. Thus, 
using (4.9a) and the fact that vanishes for large z, we compute 

: lim a3^{z, k^)ai - 



lim ^{z,k'^)= lim ^{z,k' 



1 1 



1 G{k)\ 
1 -Ij \0 F{k)J ■ 

(4.17) 

This shows that C/4 = aiA{k)ai and proves (4.9d). 

We now use continuity of the matrices $ and at the points z = ziziri to find 
U2 and U3. Let /2 denote the value of / at 2; = —iri. The conditions that ^{iC, k~^) 
be continuous at C = ri and C = are 



A 1 
fi -] 



iA{k) - U2{k)) = 



and 



h 

f2 



1 



iaiAik)ai - U3ik)) = 0, 



(4.18a) 



respectively. In view of equations (4.5), (4.6), and (4.8), the conditions that ^^^(^C) ^^) 
be continuous at C = ri and ( = — ri are 



1 



+ 2iVLh{k - ri) 



and 



-1 
1 



-1 



A{k) = 2iVlh{k-ri) 



-1 
1 



(4.18b) 



aiA{k)(Ji = 2iVth{k + n) 



-1 
1 



(4.18c) 



respectively. The top and bottom rows of each of the four matrix equations in (4.18) 
are linearly dependent since /i and /2 are purely imaginary. Combining the four 
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bottom rows into two matrix equations, we find 

Ji + 2inl{k - n) -l) ^^^^ " \2in{t- 



f2 



A{k) 

aiA{k)ai 



-1 




U2(k) 
-1 



Usik). 



^2 

Using that fi = —f2 in view of the equatorial symmetry, we deduce from these 
equations that U2{k) = Ti{k)A{k) and U^lk) = T2{k)aiA{k)ai, where Ti and T2 are 



given by (4.11). This proves (4.9b) and (4.9c). 



The properties of F and G are proved as in |LFj . 



□ 



The functions F{k) and G{k) jump across T = [—ipo, ipo]. Let F+, and , 
denote the values of F and G for k to the right and left of F, respectively. It follows 
as in |MAKj (see also [LF]) that 

^-{z,k) = ^+{z,k)D{k), /i:GF+; {z,k) = {z,k)aiD{k)ai 



where the jump matrix D is given in terms of F^ and G^ by 

D{k) 



(F+{k) 
\G+{k) 1 



F-{k) 
G-{k) 



k G F", 
(4.19) 

(4.20) 



For a given z, equation (4.19 ) provides the jump condition for a matrix RH problem 



on the Riemann surface Sz satisfied by ^{z, k). We will refer to this as the main RH 
problem]^ In general, given both the Dirichlet and Neumann boundary values for 
a BVP for the Ernst equation, it is possible to determine the spectral functions F 
and G, compute the jump matrix D, and then obtain the Ernst potential / from 
the asymptotics of the solution of the main RH problem. However, for a well-posed 



problem only one of these boundary values is specified. In our analysis of (2.6) we will 



therefore instead use the global relation and the symmetry of the boundary conditions 
to formulate an auxiliary RH problem from which F and G can be determined. 



4.3. The global relation. The equatorial symmetry of the solution / of ( 2.6 ) implies 



that the spectral functions F{k) and G{k) satisfy an important relation. Recalling 
the axis values (4.9) of the following proposition is proved in the same way as 
proposition 4.3 in |LFj . 



Proposition 4.2. The spectral functions F{k) and G{k) defined in proposition 4-1 
satisfy 



TiA+{k)aiA-^{k)T^^ = T2(JiA^{k)aiAZ^ik)aiTr 



-1 
2 ' 



keT, 



(4.21) 



where A{k) is defined in terms of F{k) and G{k) by equation {4-10) and A± denote 
the values of A to the right and left ofT, respectively. 



Equation (4.21) is referred to as the global relation. 



4.4. Linearizable boundary conditions. In general, the global relation alone is 
not sufficient for determining F and G. However, for boundary conditions satisfying 
sufficient symmetry, so-called linearizable boundary conditions, there exist another 



algebraic relation satisfied by F and G. The boundary conditions specified in (2.6) 



turn out to be linearizable. Indeed, recalling the axis values ( |4.9[ ) of the following 
proposition is proved in the same way as proposition 5.3 in [LFj . 

complete formulation of this problem also involves specifying residue conditions at the four 
points irj*^ as well as a normalization condition. 
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Proposition 4.3. The spectral functions F{k) and G{k) satisfy the relation 



{B-^K-^aia^kB){TiA+aiAl^T{^) = -{TiA+(7iA-^T:[^){B-^ K-^oia^KB), keT, 

(4.22) 

where we use the short-hand notation B and A for 

B ■■= (/(t^O) -l) ' ^ ■= '^^^+'°' ^""^ = - i^) ^ + iA;l^e-2^o(ai - I)as. 

(4.23) 

4.5. The auxiliary Riemann-Hilbert problem. Combining the relations (4.21) 
and (4.22), we can formulate a RH problem for the 2 x 2-matrix valued function M{k) 
defined by 

/ -Gik) F{k) \ 
M{k) = A{k)<j^A-\k) = f i-G(fc)^ ^J^^ 1 , A; G C. (4.24) 



Proposition 4.4. Suppose f is a solution of the BVP \2.(^ . Let fi := f{iri) G iM 
denote the value of f at z = iri. Then the spectral functions F{k) and G{k) are given 
by 

F{k)=Mi2{k), G{k)=M22{k), k£C, 
where M is the unique solution of the following RH problem: 

• M.{k) is analytic for k £ C \ {T U {— ri, ri}). 

• Across r, M{k) satisfies the jump condition 

S{k)M-{k) = -M^{k)S{k), keT, (4.25) 

where M'^ and Ai~ denote the values of M to the right and left ofV, respec- 
tively, and S{k) is defined by 

S{k) = Tf ^S-^A-Vio-gABTacJi, k£T. (4.26) 

• Ad has the asymptotic behavior 

M{k) = ai + 0{l/k), k^oo. (4.27) 

• The entries of Ai have simple poles at k = ri and k = — ri . The associated 
residues are given by 

Ke.M(k) = U-_fi ;V Re.M(k) = lL% If), (4.28) 



where 

d+ 
dC 



e2r/«) (4 29) 



and d'^ /dC, denotes the right-sided derivative. 



Proof. We deduce from (4.21) and (4.22) that the function M defined in (4.24) 
satisfies the jump condition (4.25). The asymptotic behavior (4.27) follows from the 
properties of F and G. 



By evaluating the first symmetry in ( |4.4| at the branch point iz and taking the 
limit as z approaches the regular axis, we find (cf. equations (2.63)-(2.64) in [M AKj ) 



^^^) = IWK)' ^^^)= Re/(^C) ' ^'-'"^^ 
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The poles of F and G arise since Re f{iC) = at C = Equations (4.30) together 
with the equatorial symmetry of / yield 



Res F{k) 



1 

a 



Res G{k) 



fi 



Res F{k) 



l/i 



ResG(A;) 



fi 



ri a ri a -ri a -ri a 

where a is given by ( 4.29| ). The residue conditions (4.28) follow immediately from 
these relations. 

In the last step of the proof, we show that A^2i(fe) does not have poles at the 



possible zeros of F, despite the form of (4.24). We first extend the definition (4.26) 

of S{k) to al\k£C by 

S{k) = T^^B-^ Ani+iO, k+y^aia-M-iO, k+)BT2ai. (4.31) 



This definition is consistent with ( 4.26| ). Indeed, since / is equatorially symmetric, 

" (4.32) 



Af2(z,fe+) = K^{z,-k^ 



so that A = AQ(+iO, /:;+) = Aq(— iO, A;+) for /c G F. We claim that the matrices S 
and M. satisfy 

tT{SM) = 0, A; G C. (4.33) 



Indeed, the same type of argument used to prove proposition 4.3 shows that the 
function 

n = ^^^{p + «0, /c)o-iCT3$q(p - iO, k) 
is independent of p. Evaluation at p = using the axis values yields 

n{k+) = A-^T{^B-^K^^{+i^,k+)aia^Kn{-i^,k+)BT2aiA(Ti. 
Evaluation at p = pQ yields 

Tr7^ = 0. 

The preceding two equations imply (4.33). It follows from (4.33) that — 1 must 
vanish whenever F has a zero. □ 



The auxiliary RH problem presented in proposition 4.4 can be used to determine the 
spectral functions F and G. These spectral functions can then be used to compute 
the jump matrix and to set up the main RH problem. However, in analogy with 
linearizable EVPs for other integrable PDEs, we expect that the jump condition of 
the auxiliary RH problem can also be substituted directly into the main RH problem 
with the result that the unknown quantities in the main RH problem disappear. In 
fact, an example of this mechanism was observed by Neugebauer and Meinel in the 
case of a rigidly rotating disk. They discovered that the analogs of the main and 
auxiliary RH problems can be combined into a single scalar RH problem from which 
the Ernst potential / can be directly recovered, see |MAKj . It turns out that a similar 
approach can be adopted in the present case — in the next section, we will combine the 



main and auxiliary RH problems with respective jump conditions (4.19) and (^4.25) 



into a scalar RH problem on the Riemann surface introduced in section [2| 

5. A SCALAR RiEMANN-HlLBERT PROBLEM 

We let the scalar-valued function w{k) he defined by 



w{k) 



k £ 



and define two 2 x 2-matrix valued functions £ and Q by 

C{z,k) = <l>{z,k)ai^'^{z,k), k e S^, 



(5.1) 



(5.2) 



22 

and 



JONATAN LENELLS 



Q{z, k) = -$(z, k)A{k)-^S{k)A{k)^z, k)-^ - w{k% k E S,. 

Lemma 5.1. The functions C, Q, and w have the following properties: 
• The traces and determinants of Q and C satisfy 



(5.3) 



-1, 



detQ 



-1 



w 



trQ = 0, tr£ = 0, det£ 

Q can be alternatively written as 

Q{z, k) = $(z, k)aiA{k)-^S{k)A{k)ai^{z, k)'^ + w{k)l. 
Q and C admit the symmetries 

Q{z, k') = -0-32(2;, A;+)fT3, C{z, k~) = (Jj,L{z, /c+)(T3. 

Q has no jump across T^, whereas £ satisfies the following jump conditions 
across 

{Q + wI)C. = -C+{Q + wI), k€r+, 

{Q-vjI)C- = -c+{Q-wi), ker-. 

QC = —HQ; in particular, tr(Q£) = 0. 
Lei £22 = -^^2iQii + >C22Q2i- Then 

-22 



(5.4) 

(5.5) 

(5.6) 
itions 

(5.7) 



(5. 



w has the form 



w{k) 



w^k^ + W2k'^ + wq 
^2 _ ^2 



where W4,W2,wo are real coefficients explicitly given by 
W4 = — 



W2 



Wo 



e2t/o(j)_j]^)2' 

|/oP(^7 - n^f + fM^ - 2^fe) + mri^^^h + ^1 {^rln^ - l) 

2e2C^o(J7-Oft)2 

-2ibo {f! - 2iffnnh + fi + 2innh) 



(5.9a) 
(5.9b) 
(5.9c) 



- |/o|2 {ff - AifiriQh - ^rjnl + l) + ff + ff + Aifin^h - Ar^nf^ , 
where /o and fi denote the values of f at z = +iO and z = iri, respectively. 



Proof. The first three properties in (5.4) are immediate from equations (5.1)-(5.3). 
The fourth property follows since 



det(5) 



-1 



and 



det Q = det(-$^-i(cS + wl)A<^-^) = det(5 + wl) = det5 



w 



Using (4.33) and the definitions of A and A^, a computation shows that 
A~^{S + wl)A = -aiA-^{S + wl)Aai. 



Using this identity in the definition (5.3) of Q, we find (5.5). 

The symmetries (5.6) follow from (5.2), (5.3), and (5.5) together with the first 
symmetry in (4.4). 
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By (4.19) and (4.20), ^A~^ and ^aiA~^ do not jump across F"*" and F", respec- 
tively. It follows from the expressions (5.3) and (5.5) that Q does not jump across 
F^. For k S F^, the definitions (5.2) and (5.3) of £ and Q show that 

-<!>_AZ^SM-A_<^Z^ 



-^_AZ^SA_ai^Z^ 



and 



.-1 



<^>+A^^M+SA+^^^. 



(5.10) 
(5.11) 



Using that ^A^^ does not jump across F"*" together with the jump condition (4.25), 
we see that the right-hand sides of (5.10) and (5.11) are eq ual. Similarly, using (5.5), 
the fact that ^aiA"^ has no jump across F^, and (4.25), we find the jump across 
F~. This proves (5.7). 

Since M is tracefree and tr(5A^) = 0, we deduce that MS + SM + 2wM = 0. In 
view of the definitions of Q and C, this implies QC = —HQ. 

Equation (5.8) follows by direct computation using the identity tr(Q£) = and 
the four properties in (5.4). 

The last statement concerning the form of w follows from (4.31 ) and (5.1 ) by direct 
computation. □ 



The condition that M does not jump at the endpoints of F implies that tr 5(ibz/9o) = 
0, i.e. 

two = f?Q{w2 - wapI). (5.12) 
In particular, the function h(k) defined in (2.12) vanishes at the endpoints of F. 

Lemma 5.2. There exist points {mj}f C C such that 



Q2i{k) 



(5.13) 

Proof. By (5.6), Q21 is a unique function of k, i.e. Q2i(^^) = Q2iik~). Thus 
{k'^ — rf)Q2i is an entire function of A; G C. The existence of {mj}f satisfying (5.13) 
therefore follows if we can show that 



{k' - ri)Q2i{k) 
As k ^ 00, we have 
^>(z, A;+) : 



if+f){fo+fom-nh) 



-r + o(r 



k 



1 

-1 



+ 0{l/k), A{k)=I + 0{l/k), 



00. 



00. 



(5.14) 



(5.15) 



Thus, by ([5^) and (|4.31|), 

(A:2-r2)Q(z,fc+) = - 




(/o+/o)(n-nfe)2 




-1 



f-f 2/ 
2/ / - / 



The (21)-entry of this equation yields (5.14). 
Define four points {kj}'^^-^ C C by 



^1-4^^1 + 0(^^3) 

(5.16) 

), k ^ 00. 

□ 



2 , ^in,=i(A^-fc.)(A;-fc.) 

w + 1 = — 



(A;2 



.2)2 
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We assume that the /cj's are ordered as in subsection 2.3 



Let Sz denote the double 
cover of the Riemann surface Sz defined by adding cuts j = 1, ... ,4, both 

on the upper and lower sheets of Sz- Thus a point (fc,±A, ±//) of Sz is specified by 
giving a point k G C together with a choice of sign of A and of 



1^ 



\ 



l[(k-k,){k-k,). 



We specify the sheets so that A — )■ 1 (A — — 1) as A; — )■ oo on sheets 1 and 2 (sheets 
3 and 4), and ~ /c^ (/i ~ —k^) as A; — t- oo on sheets 1 and 3 (sheets 2 and 4). As 
k crosses the cut [—iz,iz], A changes sign whereas the sign of fj, remains unchanged. 
As k crosses any of the other cuts, fi changes sign whereas the sign of A remains 
unchanged. 

Consider the function H defined by 

£22 + + 1 



where £22 = £2iQii + £22221- We fix the sign of the root V^f^ + 1 in (5.17) by 



requiring that \/w^~+T = —w^k"^ + 0{k) as — >• 00"*". Since > 0, thi s implies 
that ^/w'^ + 1 > for k S . The eigenfunction $ of the Lax pair (4.1) satisfies 
det$ = —2e'^^F{k) (see equation (2.65) in |MAK| ) . so that the entries of £ may 
have poles at the points which project to the zeros of F{k) in the Riemann fc-sphere. 



However, H has no singularities at these points. Therefore, in view of (5.8), the 
possible zeros and poles of H belong to the set in Sz which projects to {ibri}u{mj}|, 
and if H has a double pole at {mj,X,ij,), then it has a double zero at {rtij, X, —fi), 
j = l,...,4. 



By the symmetries (5.6), we have 

^22{.k, A, = C22{k, -A, n), C2i{k, A, fi) = -£2i(A;, -A, n). 

Therefore 

i/(fc,A,/x) 



1 



H{k,-X,fi) 

Similarly, we have H{k,X,fj,) = 1/H{k, X,—fi). Consequently, H — )• 1/H whenever k 



crosses one of the cuts of the two-sheeted Riemann surface defined by (2.8). We 
can therefore view H as a single- valued function on with the values on the upper 
sheet given by the values of H on sheet 1 of 5^, and the values on the lower sheet 
given by the inverses of these values. 

5.1. Formulation of the scalar RH problem. We want to formulate a scalar RH 
problem in terms of the complex-valued function ^{z, k) defined by 

^(z,fc) = ^°g^(^''\ fcGS.. (5.18) 

y 

However, since log H is a multi- valued function on S^, this definition of ip needs to be 
supplemented by a choice of branches for the logarithm. We will fix a single-valued 
representative of ip on by introducing cuts which connect the zeros and poles of H. 
Across these cuts ip will jump by multiples of 27ri/y. The problem is that even though 



(5.8) implies that all zeros and poles of H lie in the cover of the set {±^1} U {mj}f, 
the exact distribution of these zeros and poles is not known. It is therefore not clear 
at this stage how to make a consistent choice of branches. 
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We address this problem by considering the hmit in which the solution / approaches 
the Kerr solution. For a solution near the Kerr solution, we can utilize the Kerr 
expressions for F and G to compute H explicitly to first order. This will give us the 
correct choice of branches in the Kerr limit and by continuity this choice extends also 
to more general solutions. 

The Ernst potential for the Kerr black hole rotating with angular velocity Vl^ and 
with a horizon stretching from— iri to iri is given by 

^ R+e-'' + R^e''-2n 

where R± are defined by 

R± = V(±ri-C)2 + p2, 
and the parameter S G (— 7r/2,0) is related to Qh by 



2ri(l - {fkeTn2 



= 'MtrJ.' ' /i'"' = i tan(5/2) 



The value at the origin is given by 

fkerr ^ COs(<5) - 1 

■'^ cos(<5) + 1 ■ 

We consider adding a slowly rotating disk to the Kerr solution. Using the Kerr values 



for /o and /i, we compute W4,W2,wo according to (5.9) with 17 ^ 1. The branch 
points {kj}f are found by solving the equation w'^ + 1 = 0. As — t- 0, A;i and tend 
to infinity, whereas k2 and ks approach finite values. The spectral functions F^'^^^{k) 
and G^'^^^{k) are given explicitly by (cf. section 2.4 in [MAKj : note that there is a 
misprint in equation (2.349) of [MAKj ) 

_ 217|(fc^-r^) + 2ir!,/f-A:-(/f-)2 
^' 2nlik^-rl) 

For definiteness, we consider the example of ri = 1/2 and 5 = —1/2. Assuming that 
z = p + iC, with p <^ 1 and C ^ 1, we compute Q and C to first order by substituting 



the axis values (4.9a) for $ together with the values of the Kerr solution into the 



right-hand sides of (5.2) and (5.3). We find that H has double poles and double zeros 



at the points in the sets 

{mY ,^^2 , —r^ ,171^ ,rf} and {mf ,1712 , —rf ,171^ ,r^}, 

respectively. As 17 — )• 0, rrii — )■ ki and 1714 — )• k^, whereas m2 and converge to 
values close (but not equal) to k2 and k^, respectively. Let {aj, bj}f be the particular 



cycles on specified in remark 2.6, As A; traverses each of these cycles, the argument 
of H changes by the amount A arg H according to Table [T} 

A choice of branches for log H consistent with the above properties is obtained by 
introducing cuts on according to Figure |5.4[ The introduced cuts run from the 
poles of H to the zeros of H. Letting (log H)~^ and (log H)~ denote the values of log H 
for z just to the right and to the left of a cut, respectively, we have (log H)^ + 47ri = 
{logH)~. An overall choice of branch is made by requiring that logH — )• 2 log / as 



k — )• oo"*", where an appropriate branch is chosen for log/, see equation (5.20) below. 
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cycle A arg H 
ai — 47r 
02 — 47r 

Table 1. The change AargH i 
cycles in the cut basis {aj,bj}^^^. 



cycle A arg H 

63 -Utt 

64 -47r 

arg if as traverses each of the 



ki k4 




Figure 5.4. The additional cuts introduced on in order to make log if a single- 
valued function. 



For this choice of branches, we have logii(A:''') 
quently, ^p{z, k) is a unique function of k, i.e. 

i'{z,k+) = 'ip(z,k~), 

We therefore view Tp{z,-) function C — )• C. 



= - log ii(fe-), k £ C. Conse- 
yfc G C. 



Proposition 5.3. Let ip he defined by (5.18) with the choice of branches for \ogH 
specified above. Then the function 'ip{z, ■) : C ^ C has the following properties: 

• ^{z, k) is analytic for /c € C \ F U ^Uj^;^[mj, kj]^ U {±^1} 

• Across T, %l){z,k) satisfies the jump condition 

2 



i)~{z,k) = i''^{z,k) + 



VT+^it^ — w 



y{z,k+) I Vi + + w 



In 



k£T. 



(5.19) 



Across the directed intervals [kj,mj], j = 1,2; [mj,kj], j = 3,4; [ri,A;3], 
[^3)^2]; o-nd [k2,—ri], '4'{z,k) satisfies the jump condition 

ij-{z,k)=i)+{z,k)+ 

y{z,k+) 

where and denote the values of ip to the right and left of the cut, 
respectively. 

As k ^ nij, j = 1, 2, ip{z, k) satisfies 
2 

'il){z.,k) = — log(A; — mj), k ^ mj, j 



y{z,m] 



1,2. 
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log{k — mj), k^rrij, 



y{z,mj) 

As k —ri, il^{z, k) satisfies 

2 

ipiz, k) = — — log(/c + ri), 

y{z,-r{) 

As k ^ ri, ip{z,k) satisfies 

o 

'tp{z,k) 



J = 3,4. 



log(A; - ri), 



A; — 7- ri. 



A; — >■ kj, j 



L, . . . , ^, 



y{z,rj) 

As k ^ kj, j = 1, . . . ,A, ■ip(z, k) satisfies 
ip{z,k) = — , 

y 

where y = y{z,k^) for k just to the left of the cut [kj,mj] for j = 1,2, and 
just to the left of the cut [mj,kj] for j = 3,4 and is analytically continued 
around the endpoint kj so that y = y{z, k^) to the right of the cut. 
As k ^ oo, 

21og(/) 



+ 0(l/r), k^oo. 



(5.20) 



Proof. We first show that tp satisfies the jump condition (5.19). Algebraic manip- 
ulation of (5.7) using the identity tr(Q£) = and the properties in (5.4) shows that 
the functions £21 and £22 = £2iQii + £22221 satisfy 

('£21+ = -2w;£22- + (1 + 2u;2)£2i_, 
[£22+ = (1 + 2w^)C22- - 2w{l + W^)C21-, 

('£21+ = 2u;£22- + (1 + 2u;2)£2i-, 
I £22+ = (1 + 2u;2)£22_ + 2w{l + W^)C21-, 



k G r+, 
k E r-. 



(5.21) 
(5.22) 



Thus 



\/l + w'^ — w 



±2 



w 



H+{k), ker^. 



Equation (5.19) follows from here since 1 + w"^ > and V'i- + w'^ ±w>0 for k eT. 
The behavior as A; — )• itri and k — )• nij follows since H has double zeros at — ri, 

, j = 1, 2, and double poles at ri, , j = 3, 4. 
In order to find the behavior of ^ as A: — )• 00, we note that, by (5.15), 

-1 



£(z,A:+) 



/ 1 
/ -1 



0-1 



/ 1 
/ -1 



+ 0{l/k), 



Thus, by (5.16), 



£22(2;, k~^ 



4(1 + /2) 02^2 
if + f){f0 + hm - f^k 



-k^ + 0{k), 



00. 



k — ^ 00. 



Since Vw^~^^ = —Wik"^ -\- 0{k) as k ^ cxd"*", we find 
H{z,k+) = f{zf + 0{l/k), 



k — )• 00. 



28 



JONATAN LENELLS 



Therefore 



y 

□ 



5.2. Solution of the scalar RH problem. The solution of the scalar RH problem 

is 



presented in proposition 5.3 
2 

'>p{z,k) 



dk' 



ni 7p y{z, k'+){k' — k) 



\n{^Jl + w{k'f -w{k')) 



(5.23) 



dk' 



lJ[k,,m,] y{z,k'+){k' - k) 



dk' 



+ 2 



+ 



+ 



dk' 



[k,,mj] y{z,k'+){k' - k) 
fees,, 



riM] J[h,k2\ J[k2,-n]J y{z,k'+){k' -k) 

where we used that 

ln(\/l + w"^ + w) = — ln(-\/ 1 + w'^ — w), k£T. 

By deforming contours, we can replace the last three integrals on the right-hand side 
of (5.23) with an integral along 7, where 7 is the contour on defined in section [2| 
We define two divisors ^ and 5Jt on S^. ^ is defined by 

4 

ii = Y,kj, (5.24) 

whereas 371 is defined as the sum of the points in which lie above the set {mj}f 
and which are double poles of H, i.e. 



We can then write (5.23|) as 

i){z,k) 



5[R = + 1712 + + 
dk' 



(5.25) 



iri Jy+ y{z, k'){k' — k) 
dk' 



In 



(^■^l + w{k'f -w{k') 
dk' 

y{z,k'){k' -k)' 



+ 2 



^ y{z,k'){k'-k) 

where the integrals are contour integrals on and the prime on the integral from ^ 
to on indicates that the paths of integration do not necessary lie in the complement 
of the cut basis {aj, bj}. In view of (5.20), this leads to 

„/( 

log / = 



'™ k^dk 



y 



f m— - 1 

Jr+ y J-y 



k^-^dk 



y 



h{k)^^'^^ 
iT+ y 

^,,,k"'-^dk 
h{k) + 

y 



k'^dk 



7 y 

k'^-^dk 

y 



n= 1,...,4, 



(5.26a) 
(5.26b) 



where h{k) is defined by (2.12). 



Remark 5.5. Although the equations (5.26) were derived under the assumption 
that the solution is a small perturbation of the Kerr solution, they are valid more 
generally. Indeed, the crucial facts used in the derivation were that H has a double 
pole at rf and a double zero at —rf, and these properties are preserved under a 
continuous deformation. It is conceivable that the double poles of H that make up 
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OJt will change sheets under such a deformation so that (5.25) has to be modified, but 



the resulting equations (5.26) remain unchanged. 



6. ThETA FUNCTIONS 

In this section we derive explicit expressions for the Ernst potential / and the 
metric functions e^^ and a in terms of theta functions. 

6.1. Explicit expression for the Ernst potential. We will show that the right- 



hand side of (5.26a) can be expressed in terms of theta function on S^. We will 



first assume that the integration paths from to 9Jt in (5.26) lie in the fundamental 



polygon determined by the cut basis; later we will see that the result is the same also 
when this is not the case. 

Let denote the noncanonical basis of holomorphic one-forms on defined 



in (3.1) and let A be the ma trix d efined in (3.2). T hen the canonical basis is a; = A(^. 
Let u,I gC^ be defined by (|2.13l). Applying A to (|5.26bl), we find 



UJ. 



(6.1) 



Using that 



for some vector 7 G C , equation (5.26a) yields 



/ 



J. a 



-+I 



(6.2) 



where the terms involving 7 cancelled because of (5.26b). Formula (2.14) for / will 
follow if we can prove that 



e 



oc~ 00 



(6.3) 



where 0(f) := Q{v\B). Let e^'-^^ and vr^-'^ denote the jth columns of the 4x4 identity 
matrix / and the period matrix respectively. Then 

@{v + e^^^\B) = Q{v\B), Q{v + tt^^^\B) = e-2"(''^+5^^-^)G(f |S), v £ 

(6.4) 

The Jacobian Jac{T,z) of is defined as the complex torus C^/L, where L is the 
discrete lattice generated by the e*^-'''''s and the vr^-'^'s. We define the map (/? : — >• 
by 



with the contour fixed to lie within the fundamental polygon. Then ip composed 
with the projection — )• Jac(S^) is the Abel map with base point —iz. We write 
m = X;^=i Mj, where, for j = 1, . . . , 4, Mj = m+ or Mj = mj. Let /C = ip{ii) £ C^. 
An argument following pp. 322-325 in |FK| shows that /C projects to the vector of 
Riemann constants in Jac{T,z)- Thus the functions 



and 



G(^(p) - ^(mt) + /c) ' 



V^4 rP 



Pes, 



PES, 



(6.5) 



(6.6) 
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both have simple poles at the points of DJl and simple zeros at the points of K. 
Moreover, the general identity ([FK], p. 67) 



UJRS = 27Ti 



R 



j = l,...,4, 



(6.7) 



implies that as Oj is traversed the functions (6.5) and (6.6) both get multiplied by 1 



and as bj is traversed they both get multiplied by e '^'^'^ja Hence their quotient 
is a constant and we deduce that 



e((/?(p) - (^(j^) + /c)e((/7(oo-) - (^(9Ji) + fc) 
e{ip{p) - ip{m) + /c)e(^(oo-) - ip{si) + /c) ' 



Using the identity 



/ ^oo+oo- = / ^RS, i?, S G 

JS Joo- 



(6.5 



(6.9) 



and the fact that /C = </?(.ft), evaluation of (|6.8[) at P = 00''" yields 



-rf".^ + - 0(v'(oo+))e (99(00 )-/®^a;) 



0Moo+)-Xf a;)G(v9(oo-))' 
Our choice of the cut basis {aj, bj} implies that 

(P{qo^) = —ip{QG^) modulo a-periods. 



(6.10) 



(6.11) 



Thus, since Q{v) is an even function, we arrive at (6.3). 



Now suppose the integration paths from to 9Jt in (5.26) do not lie within the 
fundamental polygon. Then there exist integer vectors p,q £ such that equations 



(6.1) and (6.2) get replaced by 



u 



9Jt 



u; + Bp + q 



and 



bj '^oo■^~oo■ 



respectively. However, a computation using (6.4), (6.7), and (6.10) shows that the 



terms involving p and q cancel, so that / is still given by (2.14) 



We can now complete most of the proof of theorem 2.5; the derivation of the 
formula for e^'^ will be postponed to the appendix. We first establish the formulas 



in (2.15): The expression for Q,fi follows from (4.8); the expression for Q, follows by 



solving (5.9a) for $7 recalling that e^^" < and W4 > 0; and the expression for 
g2{/n(+j0) follows by evaluating (|2.4a|) at 2; = 



6.2. The metric functions e^^ and a. Using formula (2.14), which was established 
in the previous subsection, the expression for the metric function e^^ in (2.17) can 



be derived as follows cf. |KKSj . Since the entries of u are p urely imaginary, / G 
and e{v) = e{v) for v G (|KR2], p. 203), equation ^2.U\ yields 



G(. - n 0.) , e(^ - /; 



+ 



(6.12) 
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Let E{P,Q) denote the prime form on S^. Applying Fay's identity ( |KR2] . p. 205) 



EiP3,Pi)E{P2,Pi) 



e{v + 



P2 



e{v + 



Pi 



UJ 



Pi 



E{P^,P2)E{P^,P^) 

E{P3,Pi)E{Pi,P2) 



Pa 



Pi 



Pi 



UJ 



P2 



Q{v)e {v + 



Ps 



P2 



Pi 



UJ 



Pi 



V e 



with {Pi,P2,P3,Pi) = {iz,00- 



2Q(0)eU + 



where 



uj]Qiv + 



IZ, oo 
UJ 



to ( |6.12[ ), we find 

■•—iz 



uj\Q[v + 



UJ 



Q(o) 



= Q{v)q{v + j 
1 E{—iz, iz)E{oo~^ ,oo~) 



UJ + 



(6.13) 



UJ 



2 E{—iz, oo )E{oo~^,iz 
and we used that (lemma 3.12 in |KR2j ) 

E{—iz, oo~^)E{iz, oo~) 



E{—iz, CO )E{iz, oo^ 



1. 



By proposition 3.11 in |KR2] . Q(0) can be written as in (2.18). Letting v = u 

J^-^ UJ = u — f^-^ UJ and dividing by Q{u + f^-^ )@{u + uj), equation (6.13) 
yields 



2Q(0)- 



eiu)e{u + jrj^uj) 



UJ 



e(n + r 



e(n + n;^)G(n + /--a;) Q{u + f--uj) @{u + uj) 

Equations ( 6.12[) an d (6.14) lead to the expression for e^^ in (2.17). 
By (5.7) in |KRj, we have 



(6.14) 



(a - ao)e 



2U 



yQ{0)Q{u)Q{u + r^^^uj) 



(6.15) 



where ao € M is a constant determined by the condition that a = on the regular axis. 
In view of the expression for e^^, this yields the expression for the metric function 
a given in (2.17a). An alternative derivation of equation (6.15) is presented in the 
appendix. 

7. Axis and horizon values 



In this section we consider the limits of the formulas in theorem 2.5 as z approaches 
a point on the ^-axis. As p ^ 0, the Riemann surface degenerates since the branch 
cut [—iz,iz] shrinks to a point. This type of degeneration of a Riemann surface is 
analyzed in Chapter III of [Faj . In order to utilize the results of [Faj . we introduce 
an axis-adapted cut basis {aj,^j}f=i on by 





A 



where A 



/l -l\ 

10-1 

1-1 

\0 -I J 




The axis- adapted cut basis is displayed in Figure 7.1 in the case when (" > ReA:4. 
Note that 0,4 surrounds the coUapsing cut [—iz,iz]. According to the transformation 
formula for theta functions (Eq. (12) in [Faj), there exists a constant cq independent 
of V and B such that 

e{v\B) = coQ{v\B) (7.1) 

whenever v = A~^v = Av. We define u and I as the analogs of u and I in the 
axis-adapted basis, i.e. 



u = I hCj + oj, I = 
'r+ J-r 

Since Cj = Aco, we have u = Au. Since 



+ 



2TTiuj4, 



we have 1 = 1 — 27ritt4. Thus, introducing the shorthand notation Q(v) := Q{v\B) 



the Ernst potential (2.14) can be expressed in terms of the axis-adapted basis as 

e(^ - /_! 



-2-Kiu4+I 



e(n + n>) 



(7.2) 



Let S' denote the degenerated Riemann surface defined in (2.19). According to 
Fa] , we have the following expansions as p \,0: 

1 , 



iOj = ujj + U{p-), J = L,-z,6; i04^ —uj-^+^^ +0{p'^) 



co'j + 0{p^), J = 1,2,3; 
Bij^B'ij+0{p^), i,j = 1,2,3; 



B. 
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^Mp+^ + 0{p% 

TTl TTl 



B. 



iA 



J. + Oip" 



(7.3a) 
1,2,3, (7.3b) 
(7.3c) 



where M' G C is a constant. The path of integration from to in (7.3b) must 



be chosen as the limit of the cycle 64. For example, for Re A:3 < C < Re A;4, this path 
is [C,A;4]~ U [A;4,(^]^. Whether this path lies within the fundamental polygon on S' 
depends on the particular representatives of the homology cycles a'^ and h'^ 

7.1. Values near the regular axis. We first consider the case when z approaches 
a point on the regular axis. We define c E and c' G by 



poo roc 
/ ^, c' = / 



J. 



(7.4) 
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Lemma 7.2. The following limits hold as z approaches a point on the regular axis 
(i.e. as p — )• with ri < (): 



u 



27ri - 



2m 



(7.5a) 



e u 



@' u 



+0{p'), (7.5b) 



J —iz 



@(iL + 



CO 



e'iu' + J CO 

roo~ 

UJ 



+ 0(p2), (7.5c) 
+ 0(p2), (7.5d) 



Qi u+ I ^ + I ^ 

z J iz 



e'(n') + /3/9 + 0(p2), e{u) = Q'{u')- Pp + 0{p^), (7.5e) 



P 



e'(n' + 2c')e-^'-2^'-*^' + 0(l), 



(7.5f) 



where P £ C is a constant and u' , K' , J' , and M' are defined in (2.20), (2.21), and 
(2.29). The equations obtained from (7.5) by replacing u, u' , and J' by everywhere 
are also valid. 



Proof. The expansions in (7.5a) follow immediately from (7.3a). We will also show 
(7.5b); the proofs of the other expansions are similar. We first assume that C > Re A;4. 
In order for the argument of the theta function to have a finite limit, we shift the 
integration limit from —iz to A:4. Our choice of the cut system {dj, bj} implies that 



J —1 



Cj = c + Br + s, 



(7.6) 



where r, s G M"'' are defined by 



r = (0,0,0,1/2)^, s = (0,0,0,-1/2 



(7.7) 



Therefore, 



G u 



@{u — c — Br — s)= e^'-^a 



2TTiUN'^BN+NT{u-c~Br-s)) 



(7.8) 



Using (7.7), we can write the right-hand side of (7.8) as 



In view of (7.3c), only the terms with A'4 = and A'4 = 1 give nonzero contributions 
in the limit p — )• 0. Equations (7.3) and (7.5a) imply that the subleading terms, which 
also receive contributions from the terms with = — 1 and A'4 = 2, are of 0{p^). 
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We find 



G u 



2m[\N''^B'N'-\N''r f^^ ui'+N'T(u'-c')) 



(7.9) 



iV'i 



e' u'-c' + 



c+ 



Since 



c' + 



1 /"^^ 



+ 0{p' 



(7.10) 



this proves (7.5b) in the case when Q > ReA;4. 



Similar arguments apply when < Re/cs or Re/cs < < Re/c4. For example, if 



Refca < C < ReA;4, then equation (7.6) gets replaced by 



Cj = c + Br + t, 



where 

1 1 1 r 

Letting t' = (—1/2, —1/2, —1/2), this leads to the following analog of equation (|7.9|): 



r = (0,0, 0,1/2)^, t 



(7.11) 



coo \ / 1 /.C+ 

e( n- / w ) = e' U'-c'- 



uj'-t' 



e' u'-c' + 



\j'^\'-t^e''-^' + 0{p% (7.12) 



Taking into account that, for Re/cs < C < ReA;4, 

1 r^^ 



1 f<^ 

2 J(;- J( 



iO' , - - / Ud +t 

2 ./C- 



we again arrive at (7.5b). 



(7.13) 
□ 



Lemma 7.3. Let Q be given by (2.18). /? — )• with ri < (, 
1 



Q{u) 



e'( u' + r i^'^'' - Q' + J 



oo \ 2 



(7.14) 

The behavior of Q{0) as p ^ is given by the expression obtained by replacing u' and 
J' with zero in the right-hand side of (7.14). 



Proof. In view of ([71]), the expression for Q{u) is invariant under the change of cut 
basis from{oj,6j} to {dj,bj}, i.e. 



Q{u) 



e(n + n;^)9(n + /-'^) 



Utilizing the limits of lemma 7.2, we find (7.14). 



(7.15) 
□ 
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By applying the results of lemma 7.2 to formula ( 7.2 ) an d using that = e^"^'^*"-* = 
e^'~'^' + 0{ p^) as p — )■ 0, we find that / is given by (2.22) near the regular axis. The 
exp ression (2.24) for e^^ on the regular axis follows by applying the results of lemma 
7.3| to the equation e^*^ = ^|^e^. The limiting behavior a = 0{p'^) follows from 



]6.15); the fact that the terms of 0{p) vanish in the expansion of a is most easily 
seen from (2.2). The behavior e^*^ = 1 + 0{p^) near the regular axis follows from 
(2.17b) and the condition that k = on the regular axis; the fact that the terms of 
0{p) vanish in the expansion of e^*^ is most easily seen from (3.7). This completes 



the proof of proposition 2.8 



7.2. Values of ao and Kq. The constant oq is determined by (6.15) and the condition 
that a = on the regular axis. We find 



Oo = lim 



pe-^Q{u + !Z- + !:, 



C > n- 

Substituting into this equation the expression (|2.17al) for e^^ and passing to the 



P^o Q(0) Q{u)Q{u + PX^) 
is equa 

axis-adapted basis, we find 



p e(n + n ^ + /, 

an = iim ^, — 



oo 
iz 



OJ 



2mu4—I 



By lemmas 7.2 and 7.3, this yields 

G'(0)4 



ao 



e'(n' + 2c')^_M'-/'-2X'_ 



i\2 



Q'{u') 



(7.16) 



Expression (2.26) for ao is obtained by letting C — t- oo in (7.16). Indeed, since 

e'(o)4 



lim 



and ti', c', /' are independent of C > '"i, we find 



1 



ao = - 

The constant M' is given b"'^ 
M' = 




e i lim e-^'-^^' 



2 In X — In 4 



(7.17) 



The combination M' + 2K' remains finite in the limit C — ^ oo and we find 



lim (M' + 2K') = - hm 

(^^oo 2 i?— >oo 



R+ 



21ni?-ln4 



Equations (7.17) and (7.18) imply (2.26). 



The constant Kq in (2.17b) is determined by the condition that e 



2k 



(7.18) 



1 on the 



regular part of the axis. In order to compute Kq we first rewrite (2.17b) in terms of 
the axis-adapted cut system. Note that 



UJ 



^K+K- ~ 2vrza)4. 



(7.19) 



■'Expressions of this type are considered in |Y] 
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Since /p+ (iKi^(Ki) = 0, we find 



X e 



_ g /Ij^ , , I 'I'l 

e(o)e(/|^) 

i lim._K) if , -.ui + _— r / i(I'4-— —2 In 



C > n- 



Taking the limit as p — )• of this expression, we find (2.27). 

7.3. Values near the black hole horizon. The limits as z approaches the black 
hole horizon have a slightly different flavor than those considered in the previous 
subsection, because ua diverges as p | with < C < ''i- In fact, 

^' ^ p|0, 0<C<ri, (7.20) 



ua = -\np+ — + 0{p^ 



where 



P' 



huj'^+^^ H — lim 



7'(e) 



+ 21ne + ln4 



and 7'(e) denotes the contour 7' with the segments which lie above the interval 
[C — e, C + e] removedj^ Note that 

g2P'+2A/'^gj'^ 0<C<ri. (7.21) 

Lemma 7.4. The following limits hold as z approaches a point on the black hole 
horizon (i.e. as /> — )■ with < ( < ri): 



e u 



1 



-e'U 



^,\^J'-2M'-K' 



w + / UJ \e 



I \ 2J'~2M'~2K' 



e M+ / tD =e' U'+ / u']-Q'iu' + 



eiu+ oj]=e'iu'+ uj']+e'iu' + 



(7.22) 

+ o(i), 

a.' + y a;')e-^'+^'+0(p2), 



e U + 



-G'fn'+ / u;']e''-'''' +5 + 0ip), 



Q^u)=-Q'\^u' + uj'je'-''' +6 + Oip), 
where 5 is a constant. 



(7.23) 



Proof. We prove (7.22) in the case when Re/c3 < C, < ri] the proofs of the other 
identities are similar. For Re/cs < C < ri. 



e u 



oj = Q{u- c - Br -t), 



'''For < < ri , the contour 7' contains the covering in the upper sheet of [C, — e, ^] and the 
covering in the lower sheet of [C,, C + e]. 
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where r and t are defined in (7.11). We can write the right-hand side as 



X e 



2m 



The factor involving the divergent quantities -B44 and £14 is 

Thus the diverging terms are of 0{p^'^) and arise when A'4 = 1, 2. We find 

1 2P'-K' ST^ 2ni(^^N'TB'N'+^N'T J^;^ u'+N'T{u'~c'-t')j 



Q{ u — I uj 

' —iz 



E 



^ J_g2M'+4P'-2X' ^ ^2..(ijV'™' + |iV-/C+.'+JV-K-c'-t')) ^^^^^ 



Using (7.13) and (7.21), we find (7.22). 



□ 



Lemma 7.5. Let Q he given by {2.18 ). As p ^ G with < C < ^^i; 



Q{u) 



00 l\0 



Q(0) = '^r^ + 0{p\ 



+ 0(p^), (7.24) 



(7.25) 



Proof. By applying the limits in lemma 7.4 to (7.15), we find the statement for 



Q{u). The limit of Q{0) is obtained as in the case of the regular axis, since the 
diverging factors involving -114 are not present. □ 



By applying the limits of lemma 7.4 to equation (7.2), we find that / is given by 
(|2.30) near the horizon. Similarly, the behav ior of e'^^ near the horizon follows by 
and the expansion (7.20| to the equation e^^ 



7.5 



applying lemma 

The expression for Uhor is established by' 



Q(0) „-2TTiUi + I 



to the axis-adapted version of (6.15). 



We next show formula (2.32 



applying the results of lemmas 7.4 and 7.5 

for the behavior of e^** near the horizon. Equation 
(3.7) implies that the terms of 0{p) in the expansion of e^'^ vanish, so we only have 
to determine the leading term. Suppose < C < ^"i- T hen, recalling (7.19) and using 
that /p+ dKif (ki) = 0, we can write L as defined in (^2.16^ in terms of axis-adapted 
quantities as 



+ \ lim I / UJ + - + . . . 



21ne 1 — 2-Ki j hCj^ — 2Triui. 



In view of (|7.20|) and (|7.21|), we find 

+ 0{p), 



p^O, 0<C<ri, 



38 



JONATAN LENELLS 



of proposition 2.9 



where L' is defined by (2.28). A pplyi ng this expan sion t ogether with lemma 7.4 to 
the expression for e^'^ in theorem 2.5, we arrive at (2.32). This completes the proof 

As C 



We conclude this section by proving proposition 
whereas tends to a bounded constant. 



2.10 



we have e — t- 
The expression (2.36) for f\ := f{iri) 



therefore follows immediately from (2.23). Similarly, the statement in proposition 



2.10 regarding /o := /(+iO) follows by taking the limit p in (2.31). 



8. Parameter ranges 



In this section we consider the singularity structure of the solution (2.14) and its 
dependence on the four parameters po,ri,W2, and W4. 



8.1. Singularity structure. The solution / presented in (2.14) is continuous but 
not smooth at the point z = ziziri where the regular axis meets the horizon. Moreover, 
Im / has a jump across the disk. Away from these points / is smooth except possibly 
at points in the set where the denominator of (2.14) vanishes. Physically, we are 



interested in solutions which are singularity-free away from the disk and the horizon. 
A complete characterization of the singularity-free solutions involves determining for 
which choices of po, ri, ■W2, the set {z G P | @{u + f^-^ co) = O} is empty. We will 
not complete this analysis here, but we will indicate how a large class of singularity- 
free solutions can be constructed starting with parameters corresponding to a Kerr 
background. 



In subsection 5.1 , the Kerr solutions were parametrized in terms of the parameters 
ri > and —it/ 2 < 5 < 0. However, since the map 

2 / TT 



6 i-^ W2 



(8.1) 



tan 6 sin 6 \ 2 

is one-to-one, we may also adopt a parametrization in terms of ri > and W2 = 
2/(tan(5sin<5) > 0. Let f^^^l denote the unique Kerr solution corresponding to the 



parameters ri > and W2 > 0. Moreover, let / denote the solution in (2.14) corre- 
sponding to some strictly positive parameters po,ri,W2,W4. Then / — )• f^^^^ as po i 
and W4 ^ with r 1,11)2 held fixed. Indeed, consider perturbing a Kerr background 
solution f'^'^^^ by adding a small disk of radius pQ rotating with angular velocity In 
the limit /Jq i and O ^ 0, the jump contour in the RH problem disappears and 



Aq reduces to the identity matrix. Thus, the BVP (2.6) reduces to the Kerr black 
hole BVP and the perturbed solution / approaches f^^^ in this limit. Substituting 



the Kerr values of /o and /i and letting 17 | in (5.9), we find that 



Wo —?- 0, 



W2 



kerr\2 



1 



kerr 



tan 5 sin 5 ' 



and 



(8.2) 



as the Kerr background is approached. This leads to the relation (8.1 ) between 6 and 
W2- In view of (5.12), the vanishing limiting value of wq is achieved by letting po | 0. 
Thus, for small values of > and pn > 0, the solution / corresponding to 

fkerr t„ 
1W2 ■ 



for small values of ?i;4 > and po > 
{po,ri,W2,W4} is a small perturbation of the Kerr background solution f:^^^ 
particular, / is singularity-free for sufficiently small perturbations. By increasing W4 
and Po, larger pertubarbations of the background are obtained until the construction 
eventually breaks down and the solutions become singular. In this way, a large class 
of singularity-free solutions can be constructed. Numerical data suggest that given 
strictly positive values of the parameters po, ri, and W2, there exists an interval 



[0,w, 



maxl 
4 J 



> 0, such that all solutions / corresponding to {po,ri, 11)2,104} with 
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W4 G (0, 10™"^) are free of singularities. In the following subsection, we illustrate the 
general situation by considering a typical example. 



8.2. Dependence on W4,. We let 



Po 



n 



W2 = 3, 



(8.3) 



and consider the solution / given in (2.14) corresponding to {po, ri,W2,WA} as 11)4 > 



varies. We find that the solution is free of singularities for < 1/74 < wT""^ where 



0.27051. The example presented in section |3] corresponds to taking 1^4 



1/10. The dependence on ^4 of several parameters is displayed in Figures 8.1 and 



8.2 The parameter W4, is analogous to the variable p used in [MAKj to parametrize 



the Neugebauer-Meinel disks and the Figures 8.1 and 8.2 are the analogs of Figure 



2.9 in |MAKj . To see how the solution / becomes singular as 1^4 increases beyond 
yjmax ^ we note that as < w"^""^ increases, the ergosphere of the solution / grows 
larger and larger until it eventually, in the limit Wi f w^^""^, envelops all of spacetime. 
As Wi increases beyond w'^^^ a singularity of / enters the domain D at z = +00 and 
moves inward along the positive real axis. The graph of the singular function Re / 
1/2 > 



for W4, 



is shown in Figure 8.3 



Appendix A. Condensation of branch points 



In this appendix we show that the Ernst potential (2.14) is related via a certain 
limiting procedure to the class of solutions of the Ernst equation studied in \KM\ 
IKKS| . By applying this limiting procedure to the formula for the metric function e^'^ 
given in [KKSJ, we will also establish the expression (2.17b) for e^'' and so complete 



the proof of theorem 2.5 The limiting operation involves partially degenerating a 
Riemann surface by letting branch points coalesce along the curve r+ and at the 
points ibri. The construction of new solutions of the Ernst equation through this 
type of 'condensation' of branch points along curves was first described in [KMj . 
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0.05 0.10 0.15 0.20 0.25 



W4 



-0.4 - 



-0.8 



bo 



-1.0 



Im/i 



Figure 8.2. The dependence on of the parameters bo = Im/(-)-iO) and fi = f{iri) 



for the example specified by (8.3). 




Figure 8.3. The graph of Ref for ^4 = 1/2 > wj^ 

are too small to he visible. 



The disk and the black hole 



Let Tjz be a Riemann surface of genus 5 > 4 obtained by adding g — A brancli cuts 
{[Ej-,Fj\Y~^^ to Sz. Let ^ = —iz. Tlien is defined by the equation 

4 5-4 

-{k-i){k- Ylik - ki){k - ki) \{{k- E,){k - Fi). 



y 



i=l 



i=l 



Let {cij, be the cut basis on which is the natural generahzation of the basis 

{aj, bj}\ on S^, i.e. for j = 1, . . . , 4, aj surrounds the cut [kj, kj]; for j = 5, . . . g, aj 
surrounds the cut [Ej-4^, Fj-4]] the cycle bj enters the upper sheet on the right side of 
[—iz, iz] and exits again on the right side of [kj, kj] for j = 1, . . . , 4 and on the right 
side of [£'j_4,Fj_4] for j = 5, . . . ,g. For simplicity, we will assume that ( > Refc4. 
Let Cj = {uji, . . . ,0Jg)'^ denote the canonical dual basis and let Q{w) := Q{w\B) be 
the associated theta function. Let p,q S be vectors which are indepedent of 
z and which satisfy the reality condition Bp + q G M^. The theta function with 
characteristics p,q is defined by 



G 



p 



{v\B) = e{v + Bp + q\B)e 



2m[lp'^ Bp+p'^ [v+q)) 



Then 



e 



e 



(A.l) 
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is a solution of the Ernst equation (2.3) and the corresponding metric function e'^'^ is 
given by 



e 



e'^ = Ko- 



(0)G 



P 



men 



(A.2) 



to 



where Kq G C is a constant determined by the condition that e^"^ = 1 on the regular 
axis pCMllKKS] . 
We choose 

Ei = -ri-ie, Fi = Ei, E2 = ri - ie, F2 = E2, 
where e > is a small number, and define the Riemann surface of genus 6 by 



(k -Oik- llik - h){k - ki) ^{k - Ei){k - Fi). 



i=l 



In other words, E2 is obtained from by adding two short vertical c uts c entered at 
— ri and ri, respectively. The cut basis {dj,bj}^^i is shown in Figure 



A.l 



The condensation of branch points will now proceed in two steps: In the first step, 
we let the branch points Ej+2, Fj^2, j = 1, • • • )9 — 6 condense along the curve F. 
In doing this, the Riemann surface degenerates to and the Ernst potential / 
approaches a solution / defined in terms of theta functions on S^. Intuitively / has 
a disk, but no black hole. In the second step, we let e — )• 0. T hen degenerates to 
Sz and we will find that / approaches the solution / in (2.14). 



A.l. The first degeneration. Let 

p= {p,m) eC^ xRS-^, q = 0, 

where the components of the vector m E M^~^ satisfy < mj < 1/2, j = 1, . . . , g — 6. 
We consider the limit Ej^2,Fj+2 — ^ t^j, j = 1, • • • ,5 — 6, in which the branch cut 
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[Ejj^2, Fj^2] shrinks to a point kj G F. In this limit, (cf. Eq. (7.3)) 



{UJl, ...UJg)^ [UJl,. ..UJQ, 



1 



1 



-OJ 



2-Ki '^i ''I ' ■ ■ ■ ' 27ri %-6%-6/ ' 



Bij ^ Bij, i,j = l,...,6; 



B. 



i,j+6 



Wi, z = l,...,6, i = 1, . . . ,5 - 6; 



1 

1 



-^i+6j+6 = — In l-E'i+2 — -^1+2 1 +0(1 



j = 1, • • • ,5 - 6, i/j; 



For two points P, Q G , we have 

e 



(A.4) 



Letting N = {N,n) G x ^ and using that p = {p,m), we find that the factor 
in the sum on the right-hand side involving the diverging quantities Bj^Qj^g, j = 
l,...,g-6, is 



Consequently, since < mj < 1/2 by assumption, all terms in the sum in (A.4) 
approach zero except the ones with n = 0. We infer that the sum on the right-hand 
side of (A.4) converges to 



e 



5-6 



uj + Bp 



(A.5) 



We let the Kj's condense onto the curve F with a density determined by the measure 
dm{K) defined by 

1 dh 

dm{K) = —- — {K)dK, K G F, 

2 dK 



where h is the function defined in (|2.12|). Then, integrating by parts and using that 

00 — >• / dm{K) I Co = u, (A. 6) 



h vanishes at the endpoints of F, we find 

9-6 



where G is defined by 



u 



/ hCo. 



Combining (A.4)-(A.6), we find 



e 



CO 



e 



u 



UJ e 



(A.7) 



where L is defined by 



L 



dh 



and the prime on the integral along F indicates that the integration contour should 
be deformed slightly before evaluation so that the pole at K2 = ki is avoidedj^ 

^The result is indepedent of whether the contour is deformed to the right or to the left of the pole. 
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Applying this formula to (A.l), we arrive at the following limit of /: 

e 



f^f = 

where / G M is defined by 



e 



(A. 



r+ 



Moreover, applying equation (A.7) to the expression for e^*^ in (A.2), we find 

[^Ljrdm{^)J^^u^+^.^ (A.9) 



e(o)G(//ci.) 

where Kq is a constant independent of z. For some constant C, we have 



because both sides have simple poles at simple zeros at and are analytic elsewhere 
on E^. Hence, 



Wgf = log ' 



G TTZ + 2txiL. 



P=K- 



It follows that the last exponential factor in (A.9) is independent of z and can be 
absorbed into Kq. Thus, 



. _ G(^)e(n + /|d;) 



K, 



(A.IO) 



UJ 



A.2. The second degeneration. We now consider the degeneration of T,z as the 
cuts centered at itri collapse. In the limit e — )• 0, 

where B is the period matrix on S^, the 4x2 matrix Bi is defined by 

and the 2x2 matrix B2 is given by 

i(lne + c_) ^Ko; + - 

where c+,c_ G C are constants. Moreover, 



5-7 



2^ /r+ ^^^-^.j^ 
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e 



U + I U! 



(A.ll) 



Letting N = [N, n) G x I? and choosing, for some < a < 1/2, 

p= (0,0, 0,0, a, -a), 



the same type of argument that led to (A. 5 ) shows that ah terms in the sum approach 



as e — )• except those with n = 0. It fohows that the sum in ( A.ll ) converges to 

L 



( n + / a; I where u 



hoj + a 













On the other hand, the rightmost exponential factor in (A.ll) can be written as 



^2'Ki{lpTBp+pT(u+J^ Co)) ^ g7ria2 (^55 -2556+^66) 



X e 



Applying these formulas to (A.8), we find 



(A.12) 



where 



+ a 



r+ 



The contours in the integrals J 1 and f }_ are the limits of the cycles 65 and Lq, 
respectively. By deforming these contours, we find that 



00^00 5 



UJ = I u, 

1 / --7 



where 7 is the contour on defined in section [2] Therefore, in the limit a — t- 1/2, 
u and / become exactly the u and / of theorem |2.5| and the solution / in (A.12) 



becomes the Ernst potential in (2.14). This provides the promised link between the 



solutions in |KM| and the solution presented in this paper. 



By applying the same limiting procedure to equation (A. 10) we will determine the 
corresponding metric function e'^'^. We have 



0. 



Hence, for the quotient in (A. 10) we find 

Kq — = = = e e ''ill 



e(o)e(//d)) 



0(o)e(//' 



(A.13) 



X e 



1 "^1 '^I'^i e 1 1 1 1 (1 + C(e)) 



This expression vanishes in the limit e — )• 0. However, this behavior is compensated 
by the fact that the constant Kq diverges as e — )• 0, so that the limit of e^'^ is finite 



and non-zero. The last exponential factor on the right-hand side of (A.13) can be 
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absorbed into Kq. Indeed, the same type of argument that we used to find (A. 10) 



shows that this factor is independent of z. The constants c+ and c_ are given by the 



expressions obtained by replacing C, in the right-hand side of (2.29) with ri and —ri, 
respectively. Letting a — )• 1/2 and using that 

/ /■(-ri-<5)+ r{ri-S)+ \ 

/ + - -4 In (5 - 



lim 

5^0 



+ 



-ri-5)- 



2 hm 



'71 (-5) 



21n5 



72(5) 



we infer that the limit of e is given by (2.17b) 



We finally point out that formula (6.15 ) for the metric function a can be derived in 



a similar way. Indeed, the metric function a corresponding to the solution / in (A.l) 
is given by |KKS] 



(a - ao)e 



2U 



(e 


p 
1. 


(0)G 


p 




^0(0)6 


P 




p 

_q_ 





where 



Q(o) 



An application of the above limiting procedure to this expression yields (6.15) 
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